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Limiting conditions for soot particle inception in spherical diffusion flames were investigated 
numerically. The flames were modeled using a one-dimensional, time accurate diffusion 
flame code with detailed chemistry and transport and an optically thick radiation model. 
Seventeen normal and inverse flames were considered, covering a wide range of 
stoichiometric mixture fraction, adiabatic flame temperature, residence time and scalar 
dissipation rate. These flames were previously observed to reach their sooting limits after 2 s 
of microgravity. Sooting-limit diffusion flames with scalar dissipation rate lower than 2 s' 1 
were found to have temperatures near 1400 K where C/O = 0.51, whereas flames with 
greater scalar dissipation rate required increased temperatures. This finding was valid 
across a broad range of fuel and oxidizer compositions and convection directions. 


Nomenclature 


A 

area 

C/O 

carbon to oxygen atom number ratio 

c/o crit 

minimum C/O for soot formation 

c p 

constant-pressure heat capacity of the gas mixture 

Cp.lc 

constant-pressure heat capacity of species k 

D 

mixture-averaged diffusion coefficient 

h 

specific enthalpy 

k 

species number 

K 

total number of species 

m 

mass flow rate 

M 

molecular weight 

r 

radius; oxygen/fuel mass ratio, Eq. (8) 

T 

temperature 
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Tf peak temperature 

t res residence time 

V k diffusion velocity of species k 

V(r) velocity at radius r 

Y mass fraction 

Z mixture fraction 

Z st stoichiometric mixture fraction 

1 mixture thermal conductivity 

p gas density 

X scalar dissipation rate, Eq. (13) 

co molar chemical production rate 

I. Introduction 

F ormation of soot in combustion is an active research topic. 1 " 3 One fundamental measure of flame sooting 
behavior is sooting limits. Sooting limits of spherical microgravity diffusion flames were observed by 
Sunderland et al. 4 The present work is an investigation into those flames with the aid of a detailed computational 
model. 

In the past, most data on fundamental sooting limits came from studies of laminar premixed flames. 1 ' 2 ' 5 " 8 . One 
reason for this is that both temperature and carbon-to-oxygen atom ratio, C/O, are nearly constant in the soot- 
forming regions of premixed flames. Past studies found sooting limits in premixed flames to arise from a 
competition between fuel pyrolysis and oxidation of soot precursors. 8 It was found that soot inception in premixed 
flames cannot occur when C/O is below a critical value, about 0.6 for ethylene. 1 ' 2 ' 5 ’ 7 ’ 8 

A similar competition can occur in diffusion flames on the fuel side owing to the presence of oxygen in species 
such as C0 2 and H 2 0. 4,9 Du et al. 9 showed that adding C0 2 to the fuel side of diffusion flames can suppress soot 
formation chemically. Despite the differences between soot inception in premixed and nonpremixed flames, the C/O 
atom ratio has been shown to be relevant to sooting limits in diffusion flames. 4 ’ 10 Past experimental work on 
spherical diffusion flames in microgravity identified a critical local C/O value of 0.59 for ethylene. 4 

Temperature plays a different role in premixed and nonpremixed flames: increasing temperature suppresses soot 
formation in premixed flames, whereas it enhances soot formation in nonpremixed flames. 2 ' 8 Previous research on 
diffusion flames has identified an onset temperature at which soot particles are first observed to be in the range of 
1250- 1650 K. 2 ' 11 " 14 

In addition to C/O ratio and temperature, residence time (or scalar dissipation rate) is the third parameter critical 
to soot inception limits. Soot induction times of 0.8 - 15 ms were reported by Tesner and Shurupov in Ref. 15 for 
acetylene/nitrogen mixtures at 1473 K. Strain rates of 30 - 200 s" 1 were observed to prevent soot formation in 
counterflow diffusion flames. 16 ’ 17 

Microgravity allows the observation of strain-free diffusion flames. It also allows conditions with much longer 
residence times than can be obtained in normal gravity. These long residence times can yield new insights into soot 
inception limits, but they cause increased radiative losses that require consideration. 

This study seeks to further investigate the effects of local C/O atom ratio, local temperature, and residence time 
on sooting limits of spherical diffusion flames. The flames under consideration are the 17 microgravity flames of 
Ref. 4, having both normal and inverse convection directions and disparate stoichiometric mixture fraction, adiabatic 
temperature, and residence time. The model used here is a one-dimensional, time accurate diffusion flame code with 
detailed chemistry and transport and an optically thick radiation model. 

II. Numerical model 

The conservation equations were solved numerically using a flame code that includes detailed kinetics, based on 
GRI Mech 3.0 and involving 53 species and 325 reactions, and detailed transport properties. The numerical code 
was modified from the Sandia burner-stabilized PREMIX code, 18 which was originally developed for the study of 
one-dimensional freely propagating and burner-stabilized premixed laminar flames. The code was adapted for 
spherical diffusion flames allowing for optically thick radiative heat losses. 19 Both steady state and transient flames 
can be modeled. 

Modifications were introduced to improve agreement with experimental results. Species and heat diffusivity has 
been increased by 30% from the original version." 0 This increase yields higher peak temperatures and smaller 
flames. 

In transient cases, the conservation of mass, energy, and gas species equations solved by the code are: 
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(1) 



The equations were solved subject to the following boundary conditions: 

r = r b :T = T b ;Y k = ^ , k = \,2,—,K and (4) 

m b 

r — » +oo ; T — » T^Yk — » F^ )00 , k = 1,2,...,* , (5) 

where subscripts /> and oo refer to conditions at the burner surface and at the ambient, respectively. Variable T fc() is 
the mass fraction of the k th species in the burner gas supply. Because drop tower experiments did not show a 
significant increase in burner surface temperature, 19 T h was taken to be constant. The outer boundary temperature 
(295 K), species mass fractions, and pressure are held constant and specified as input data. However the velocity at 
the outer boundary is finite. 

The computation requires a set of initial conditions that resemble the experimental ignition process. Following 
the approach adopted by Tse et al., 21 the initial (ignition) conditions were prescribed as the steady-state solution of 
the same flame in a compressed domain. These ignition conditions had no radiation and assumed adiabatic 
conditions at the burner exit and a constant temperature at the outer boundary. The compressed domain was chosen 
to extend 1.2 cm from the burner center. The steady-state computations begin with a set of prescribed initial 
distributions of grid points, temperature, species concentrations, and an estimate of the location and thickness of the 
reaction region. After the steady-state solution is obtained, the grid is adapted to reduce the gradients and curvature 
and to improve accuracy. New grid points are added until all values of the gradients and curvature are below the user 
specified limits, and further addition of grid points does not affect the solution. 

The steady state solution on the compressed domain is then used as the initial condition of the transient 
computations by expanding the domain of computation to 100 cm from the center of the burner, by adding grid 
points to fill the gap between the compressed and normal domains. All the computations used a grid of about 200 
mesh points, finer in the area of high gradients (flame location), becoming coarser close to the outer boundary. Test 
cases were run to insure grid independence. Typically, about 120 extra mesh points were added from the compressed 
to the extended computation domain. The initial conditions at those extra mesh points were taken to be the user 
specified values for the outer boundary. 

For the transient computation of the flame structure, radiative losses were considered. Radiation was considered 
optically thick and caused only by the participation of C0 2 , H 2 0, and CO. The radiative properties of these gases 
were formulated by a statistical narrow-band model with a spectral bandwidth of 25 cm" 1 . The emissivities were 
extracted from the line-by-line values given by the HITRAN database. 24 


III. Experimental 

The tests considered here were performed and reported in Ref. 4. Further details of the apparatus are contained in 
Ref. 22. These experiments were conducted in microgravity in the NASA Glenn 2.2 s drop tower. The burner was a 
6.4 mm diameter porous stainless steel sphere. All tests were conducted in quiescent ambient gas at 295 K and 0.98 
bar and ignition was performed in microgravity. The tests employed three gases: ethylene, nitrogen, and oxygen. A 
color video camera was used to image the flames. 

A summary of the 17 sooting limit flames is given in Table 1. Burner flow rates were selected such that all 
flames involved a steady-state ethylene consumption rate of 1.51 mg/s and a heat release rate of 71 W, assuming 
complete combustion. Normal flames are realized when fuel discharges from the burner into an oxidizer atmosphere, 
while inverse flames are realized when oxidizer flows from the burner into an atmosphere of fuel. The fuel and 
oxygen mole fractions in the supply gases, X C 2 H 4 .o and X 0 2 ,o, vary widely. This yielded a wide range of 
stoichiometric mixture fraction, Z s ,, as shown. Adiabatic flame temperatures were calculated using Chemical 
Equilibrium with Applications, CEA. 23 
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Table 1. Summary of the sooting limit flames, reproduced from Ref. 4. 


Flame 

Environment 

Vc2H4.0 

V02.0 


^resi S 

Tad, K 

Tr, K 

1 

Oxidizer 

1 

0.22 

0.065 

2.72 

2390 

1545 

2 

Oxidizer 

0.6 

0.21 

0.102 

1.63 

2326 

1492 

3 

Oxidizer 

0.31 

0.21 

0.18 

0.91 

2226 

1479 

4 

Oxidizer 

0.25 

0.23 

0.225 

0.665 

2238 

1498 

5 

Oxidizer 

0.18 

0.28 

0.333 

0.351 

2306 

1592 

6 

Oxidizer 

0.17 

0.29 

0.353 

0.33 

2308 

1593 

7 

Oxidizer 

0.11 

0.5 

0.586 

0.11 

2381 

1795 

8 

Oxidizer 

0.11 

0.8 

0.685 

0.044 

2528 

2057 

9 

Oxidizer 

0.15 

1 

0.661 

0.024 

2740 

2262 

10 

Fuel 

1 

0.13 

0.041 

0.059 

1847 

1581 

11 

Fuel 

0.8 

0.13 

0.051 

0.072 

1835 

1549 

12 

Fuel 

0.6 

0.13 

0.066 

0.086 

1814 

1515 

13 

Fuel 

0.21 

0.25 

0.277 

0.119 

2274 

1689 

14 

Fuel 

0.19 

0.3 

0.336 

0.122 

2370 

1736 

15 

Fuel 

0.15 

0.5 

0.509 

0.148 

2539 

1802 

16 

Fuel 

0.12 

0.8 

0.666 

0.279 

2578 

1729 

17 

Fuel 

0.13 

1 

0.692 

0.249 

2670 

1814 


Four representative flames at or near their sooting limits are shown in Fig. 1. These flames represent both 
convection toward oxidizer (normal flames) and convection toward fuel (inverse flames). The flames of Figures lb 

and Id are considered here to be at experimental 


(a) 1 8% C 2 H 4 -> 27% 0 2 (b) 1 8% C 2 H 4 28% 0 2 


4 


30 mm 


OO 


(c) 0 2 -> 12%C 2 H 4 (d) 0 2 — > 13% C 2 H 4 


Figure 1. Color images of representatives flames below 
the sooting limits (a and c) and at sooting limits (b and d) 
for convection toward oxidizer (a and b) and convection 
toward fuel (c and d). Flames (b) and (d) correspond to 
Flames 5 and 17 in Table 1. Images were taken just 
before drop termination and are reproduced from Ref. 4. 


sooting limits because a small reduction in 
reactant concentration yields blue conditions at 
drop end. Note that soot, when present, appears 
inside the flame sheet for normal flames and 
outside the flame sheet for inverse flames. 

Flame temperatures 7), taken as the peak 
temperature at 2 s, are obtained from the 
numerical model. Residence time is the time for a 
parcel of gas to convect from the burner surface to 
the flame sheet. It is obtained by integrating the 
predicted flow field at 2 s: 

flame 


trpv 


- 1 


dr 

v{7) 


(6) 


burner 

These residence times (see Table 1) are slightly 
different from the more approximate values 
obtained in Ref. 4 from a simple scaling analysis. 

IV. Results 

The 17 spherical diffusion flames summarized in 
Table 1 were simulated using the modified 
PREMIX code. Predictions at 2 s after ignition 
were emphasized to better understand the 
conditions associated with the observed sooting 
limits. The flames modeled covered both normal 
and inverse configurations for a wide range of 
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Figure 2. Predicted flame radius and peak temperature of Flame 10, 
from ignition (t=0 s) to 3 s. 


residence time, Z st and adiabatic flame 
temperature. 

Figure 2 presents the Flame 10 
computed flame radius and peak 
temperature evolution. Flame radius is 
defined as the peak temperature 
radius. 

The proximity of the cold outer 
boundary in the compressed (ignition) 
domain promotes conductive heat 
losses, reducing the peak temperature 
below its adiabatic value. 
Immediately after ignition, the 
absence of conductive heat losses at 
the outer boundary increases slightly 
the flame temperature. The radius 
decreases slightly before increasing 
again. 

The measurements and 
computations indicate that neither 
flame size nor peak temperature 
reaches steady state within 2 s for any 
of the present flames. Flame 


expansion causes temperatures to drop owing to increased radiative losses, which are proportional to flame surface 
area 19 . At 2 s, the radiative heat loss fractions are around 0.42 for all flames except Flames 10, 11 and 12, whose 
heat loss fractions are about 0.28. Flame expansion is a critical aspect of all these flames. 

A similar evolution occurs for flame temperature. For the flames of Table 1, peak flame temperatures at 2 s span 
a broad range of 1479 - 2262 K. These temperatures are too high to cause radiative extinction at 2 s. 19 

The calculation of mixture fraction, Z, for these flames is complicated by the presence of N 2 in both the fuel and 
oxidizer in most cases. Furthermore, even for cases without fuel dilution, the conventional definition of 
Z CH = Y c + Y h does not yield peak temperatures at Z st . As discussed below, these complications can be overcome by 
including Y 0 in the defined mixture fraction. 25 

A definition of Z, based on the mass fractions of C, O, and H atoms, was derived by Bilger 25 for methane. For 
any hydrocarbon fuel, C„F1„„ this becomes 


z CHO = 


Y C n H m ~ Y Cn H m .2 + H y 02 .2 " Y 0 2 ) 


Y C n H m .1 Y(2 jj 2 + 


-(y 02 ,2-Yo 2 , l)’ 


where 1 and 2 denote fuel and oxidizer supply, respectively, and where 

( n + m/4 )M 0l 

r = — . 

M C„H m 

One can express Y CmHn as a function of Y c and Y H .This gives 

nM 


(7) 


( 8 ) 


Y C =- 


Yh = 


m c„h„ 
nM 
Mr 


-Yr H 
v n n m 


= 0.5 Mr* u 
'-n n m 


-Y C // 


l c 


y h 

nM fj mM(j 


(9) 


l C n H m 

For 0 2 , we have the following relation: 


Inserting Eqs. (8 - 10) into Eq. (7) yields 


*02 = 


M Q 2 
2 M 0 


r O 


(10) 
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(11) 


Z CHO = 


Yh~Y H ,2 ] Y c - y c,2 ] Yq2-Yq 

mM fj nM(j ( n + m/4)MQ 

y H\ ~ y H, 2 + Y C \ ~ Y c 2 | Yq 2 ~ Yq i 
mM uMq (n + m/4)M(j 


which, for ethylene, becomes: 


z CHO = 


Y H ~ Y H ,2 | Y c ~ Yq 2 | Yq 2 ~ Yq 
4 M h 2 M c 3 Mq 

Y H\ ~ Y H, 2 | Y C \ ~ Y Ci2 | Yq 2 ~ Yq i 
4 M h 2 M c 2 >Mq 


( 12 ) 


Zhou and Mahalingam 26 evaluated various definitions of Z by plotting predicted temperature profiles with 
respect to mixture fraction. This approach is shown in Fig. 3 for Flame 10. Two definitions ofZ are considered, each 
yielding peak temperatures at different values of Z. The Z CH definition results in peak temperatures on the rich side 
of the known Z st of 0.041. In agreement with Ref. 26, the Z C ho definition of Eq. (12) does not have this problem; it 
therefore is adopted for the results that follow. 

Figure 4 depicts the predicted structure of Flame 10 at 2 s after ignition. The inverse configuration of this flame 
is evident in the increase of C/O with increasing radius. The temperature field extends to a radius of 9 cm at time 2 s. 
All 17 flames present higher temperature and species concentration gradients inside the flames than outside. 

Past work 4 proposed that with sufficiently long residence times, at its sooting limit each flame has a location 
where C/O and T reach critical values below which soot cannot form. The numerical model allows this hypothesis to 
be tested using the present flames, which were observed to reach sooting limits at 2 s. Figure 4 reveals that for 
Flame 10 there is a location on the fuel side where C/O is 0.51 and T is about 1400 K. (The selection of these critical 
values is explained below.) Similar behavior with respect to these critical values is observed in all the present flames 
that have sufficiently low scalar dissipation rates, as shown below. 

Although spherical diffusion flames are strain-free, their scalar dissipation rates are nonzero. As with strain, 
scalar dissipation rate can dramatically affect soot inception. Scalar dissipation rate is defined here as 


Z = 2 D N2 ( dZ / d/' ) 2 . (13) 

The diffusivity coefficient used in Eq. (13) is the diffusivity of N 2 into the multicomponent mixture at the local 
temperature. The choice was made because N 2 concentrations are relatively high near the peak temperature locations 
for all 17 flames. 

Temperature and C/O ratio profiles, similar to those of Fig. 4, were prepared for each of the present flames. 



Radius (cm) 



10 


1 


0.1 


0.01 


0.001 


Figure 3. Predicted temperatures of Flame 10 
plotted with respect to mixture fraction. Two 
definitions of mixture fraction are shown: Z c h and 
Zcho- Results are at 2 s after ignition. 


Figure 4. Predicted radial profiles of temperature, 
mixture fraction, scalar dissipation rate, and C/O 
ratio for Flame 10. The burner surface is at 
0.3175 cm. Results are at 2 s after ignition. 


6 

American Institute of Aeronautics and Astronautics 


s) /- ‘0/3 ‘z 





Various C/O ratios near 0.6 were then considered. For each of these, the associated local temperature, T(C/0), and 
local scalar dissipation rate, /(C/O) were found for each flame. For each value of C/O, the standard deviation of 
these temperatures was considered. Flames 7-9 had the highest /(C/O) and thus were excluded from the standard 
deviation calculations. 

Figure 5 shows the relationship between the standard deviation of 7(C/0) and C/O. The minimum occurs at a 
C/O of 0.51. The standard deviation is not strongly dependent on C/O for C/O in the range of 0.45 - 0.6. Based on 
these results, the present flames indicate a critical C/O ratio of C/O cr ; t = 0.51. This is in reasonable agreement with 
previous reports of a value of 0.6 in ethylene premixed flames 1 ’ 2 ’ 5 ' 7 ’ 8 , a value of 0.59 obtained in a simpler analysis of 
these spherical flames, 4 and a value of 0.6 identified for gas-jet diffusion flames. 27 

Predicated on a value of C/O cr i t = 0.51, the associated local critical temperature, T crih and local scalar dissipation 
rate for each flame are shown in Figure 6. The ordinate here is inverse scalar dissipation rate, which indicates a 
characteristic mixing time. The T crit values are similar for most of the present flames. The exceptions are Flames 7 - 
9, which have the highest local scalar dissipation rates, i.e., l//(C/O cr j t ) < 0.5 s, and thus require higher local 

temperatures to reach their sooting limits. Excluding 
Flames 7-9, the average T crit is 1400 K. 

The inverse of scalar dissipation rate gives an 
estimate of the characteristic mixing time. At the soot 
inception location, this time has a broad range, from 
0.13 s (Flame 9) to 15.6 s (Flame 1). This characteristic 
time scales roughly with residence time for the normal 
flames and with inverse residence time for the inverse 
flames. Thus residence time is not a good indicator of 
scalar dissipation rate in spherical inverse flames. 

Figure 7 shows the critical temperatures of Fig. 6 
plotted with respect to Z st . The critical temperature for 
soot formation is seen to be independent of structure 
(represented by Z, t ) when Flames 7 - 9 are excluded. 
These are identifiable as the three flames in Fig. 7 with 
the highest temperatures (and open symbols). 

The present 17 sooting limit flames represent a 
C/O Ratio broad range of conditions in terms of fuel and oxidizer 

Figure 5. Standard deviation of T(C/0) plotted with supply concentrations, Z st , peak temperature, adiabatic 

respect to C/O. Flames 7-9 had the highest /(C/O) flame temperature, convection direction, residence 

and thus were excluded. time, and scalar dissipation rate. Nevertheless, when 
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Figure 6. Predicted critical temperatures, at C/O crit = 
0.51, versus inverse scalar dissipation rate at this 
location for the present 17 flames. 



Figure 7. Predicted critical temperature, at C/O crit = 
0.51, versus Z s , for the present 17 flames. 
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the three flames with the highest local scalar dissipation rate are excluded, a critical C/O ratio of 0.51, and a critical 
temperature of 1400 K, occur at the sooting limits. 

V. Conclusions 

Computations with detailed chemistry and transport of spherical diffusion flames at sooting limits have been 
performed. Seventeen flames were considered, these having been previously observed to reach sooting limits after 
2 s of microgravity. The key findings are as follows: 

1) A mixture fraction calculated from local mass fractions of carbon and hydrogen atoms is not valid for these 
flames. An improved definition that also includes oxygen atom results in peak temperature near the 
stoichiometric mixture fraction. 

2) Sooting-limit diffusion flames with sufficiently low scalar dissipation rate have similar temperatures where C/O 
= 0.51. For flames with scalar dissipation rate at this location lower than 2 s' 1 , these temperatures average 
1400 K and are not affected by Z sl . 

3) Flames with scalar dissipation rate greater than 2 s" 1 at this location require increased temperatures. 
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